Spectral signatures of the Luttinger liquid to charge-density- wave transition 
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Electron- and phonon spectral functions of the one-dimensional, spinless-fermion Holstein model at half 
filling are calculated in the four distinct regimes of the phase diagram, corresponding to an attractive or repulsive 
Luttinger liquid at weak electron-phonon coupling, and a band- or polaronic insulator at strong coupling. The 
results obtained by means of kernel polynomial and systematic cluster approaches reveal substantially different 
physics in these regimes and further indicate that the size of the phonon frequency significantly affects the nature 
of the quantum Peierls phase transition, the latter being either of the soft-mode or central-peak type. The generic 
features observed are relevant to several classes of low-dimensional materials. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.30.+h, 71.45.Lr, 71.38.-k 



I 

C 

o 
o 



(N 
> 

m 

o 
o 

03 



i 

o 
o 



X 



I. INTRODUCTION 

Low-dimensional materials like halogen-bridged transition 
metal complexes, ferroelectric perovskites, conjugated poly- 
mers, or organic charge-transfer salts are very susceptible to 
structural distortions driven by electron-phonon (EP) interac- 
tion. At commensurate band fillings these systems might un- 
dergo a Peierls or charge-density-wave (CDW) instability ac- 
companied by a dimerization of the lattice but, unlike in con- 
ventional metals, both quantum lattice fluctuations and strong 
electronic correlations are important. 1 The challenge of under- 
standing the related metal-insulator transition, especially in 
the strong-EP coupling regime and in conjunction with strong 
electronic correlations, has renewed the interest in models of 
interacting electrons and phonons^ii^ 

In this work, we compare in detail the spectral signatures 
of four physically distinct regimes in the phase diagram of 
the one-dimensional (ID) Holstein model of spinless fermions 
(HMSF) (Fig.Q. To this end, we calculate the single-particle 
spectral functions of electrons and phonons, using exact di- 
agonalization (ED) in combination with cluster perturbation 
theory (CPT)& and the kernel polynomial method (KPM). 7 
Besides the influence of the degree of phonon retardation on 
the nature of the Peierls transition, the wave-number depen- 
dence of the spectral functions missed in previous dynamical 
mean-field theory studies 3 turns out to be crucial for the ID 
case considered here. 

The paper is organized as follows. The model is introduced 
and its phase diagram reviewed in Sec. In Sec. [TO] we 
present the methods used, and results are discussed in Sec. lIVI 
Finally, we summarize in Sec.lVl 



II. MODEL 
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is particularly rewarding to study a number of interesting gen- 
eral phenomena because it exhibits at half filling, as the EP 
coupling increases, a zero-temperature quantum phase tran- 
sition from a metallic Luttinger liquid (LL) to an insulating 
Peierls state. In Eq. iff), cj {b\) creates an electron (disper- 
sionless phonon of energy uq) at site i of a ID lattice with 
N sites. The first term describes hopping processes (oc t) 
between neighboring sites the second term gives the 

elastic and kinetic lattice energy, and the third term accounts 
for the local coupling (oc g) between the lattice displacement 
ii = b\ + b i and the electron density hi = c\c { , all in the limit 
of infinite on-site (Hubbard) Coulomb repulsion, i.e., the site 
occupation numbers are ni = 0,1. Important parameters to be 
used here are the adiabaticity ratio coo/t and the dimensionless 
coupling A = £p/2t, where e v — g 2 u>o is the polaron binding 
energy. 



Despite its simplicity, the HMSF is not exactly solvable 
and a wide range of analytical and numerical methods have 
been applied to map out the ground-state phase diagram in 
the g—uio plane (see references in Ref. 8). At present the 
probably most precise phase boundary is obtained by ex- 
act diagonalization and density matrix renormalization group 
(DMRG) techniques. 2 ' 9 More recent large-scale DMRG cal- 
culations supplemented by a finite-size analysis have proved 
that at low (high) phonon frequencies the metallic LL phase 
is characterised by an attractive (repulsive) interaction^ But 
also above g c (wo), where long-range CDW order sets in, there 
exist two physically distinct regimes, which can be classified, 
e.g., by their different optical response, 9 either as a band in- 
sulator in the adiabatic regime u>o/t <C 1 or as a polaronic su- 
perlattice in the limit of large phonon frequencies u>o/t 3> 1 
(see Fig.Q. Including the spin degrees of freedom and a finite 
Hubbard interaction allows for additional quantum phase tran- 
sitions between Peierls and Mott insulating phases. 4 The half- 
filled HMSF captures the relevant physics of the more gen- 
eral quarter-filled Hubbard-Holstein model 3 4 in the regime of 
large Hubbard repulsion U ~^>t, often realized in experiment, 
where on-site bipolaron formation is suppressed. 
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FIG. 1: DMRG ground-state phase diagram of the ID half-filled 
HMSF (cf. Refs. 2 9). In the adiabatic limit lo — ► the critical 
dimensionless coupling constant A c converges to zero. For luq > 0, 
due to quantum phonon fluctuations, there exists a finite critical cou- 
pling g c (luq). For luo — > oo, the model exhibits a Kosterlitz-Thouless 
phase transition near g* — g c (uJo — * oo). 2 Crosses indicate parame- 
ter sets used below. 
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The T = phonon spectral function is defined as 

B(q,Lu) = --lmD(q,Lu) 

IT 



with 



D{q,uj)= \im {tp \x q — — -X- q \ip ) 

?7^0+ LU + 17] — H 
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for lu > and x q = N' 1 / 2 £\ Xje^ri. For the HMSF Q, 
w) is symmetric in g, and we have a dispersionless bare 
propagator Do(q, uS) = 2luq/(lu 2 — u>q). EP interaction will 
renormalize the phonon frequency, whereby D(q,ui) attains 
a g-dependence. Concerning notation, we shall reserve the 
capital Q for allowed wavenumbers of finite clusters. 

IV. RESULTS 



III. METHODS 

The T — electron spectral function is related to the one- 
electron Green function via 

A(k,u) = — Im G(k,uA = A + (k,cj) + A~(k ) u>) , (2) 

7T 

where 

A ± (k,u;) = --Im — 1— - c±|^) , (3) 

7T r;^0+ U + ITj =p il 

= c fc , = ct, and |-0o) 1S tne ground state of the HMSF 
with Af e = N/2 electrons. A^(k,uj) [A + (k,uj)] describes 
(inverse) photoemission of an (injected) electron with mo- 
mentum k and energy uo. 

Relying on approximation-free numerical approaches (here 
ED) we can calculate exact Green functions for finite systems 
only. To obtain an approximation of G(k,u>) for the infinite 
lattice, we can exploit CPT 6 To this end, we divide the infi- 
nite system into identical clusters of N c sites each, and deter- 
mine the electron cluster Green function g\j (u) for all non- 
equivalent pairs of sites i, j = 1, . . . , N c by the KPM (for de- 
tails see Ref. 7). The phonon Hilbert space is truncated such 
that the resulting error of the spectra is negligible (< 10 -4 ), 
and we have used 1024 Chebyshev moments. In a second step, 
the Green function G(k, u>) for the infinite lattice is obtained 
from the first-order result 6 



Gij(k,u>) 
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(4) 



of a strong-coupling expansion in the inter-cluster hopping op- 



The crosses in the phase diagram (Fig. ^ mark the four 
parameter sets to be used in the sequel. 

Figures|2ja) and (c) show results for the CPT electron spec- 
tral functions A^ (k, lu) in the metallic region. In the adiabatic 
(attractive LL) regime [Fig.[2a)], we find a rather pronounced 
peak at the Fermi level Ep but nevertheless the system is not 
a Fermi liquid. The non-universal LL parameters (charge ve- 
locity and interaction coefficient) can be determined from a 
DMRG finite-size scaling. 10 At k = kp, very little spectral 
weight is contained in the incoherent part of the spectrum. By 
contrast, away from the Fermi momentum, almost all of the 
spectral weight resides in the incoherent part, whose maxi- 
mum follows quite closely the free dispersion —2tcosk [no- 
tice the symmetry of A~ (k < kp) and A + (k > kp)]. As these 
states are accessible only via (multi-) phonon excitations, the 
width of the incoherent band is proportional to g 2 . 

The spectrum in the anti-adiabatic (repulsive LL) regime, 
Fig.|2jc), looks significantly different. Due to the larger cou- 
pling A = e p /2t = 1 (with respect to the bare electronic 
bandwidth), as compared to A = 0.5 in Fig. Efa), the carri- 
ers are more strongly renormalized. This is reflected in the 
reduced bandwidth W w 3t, with the spectral weight being 
more evenly distributed than in the adiabatic case. The fact 
that luq > W gives rise to multiple, non-dispersive side-bands 
at energies Ep ± tiloq. 

In Figs.|2jb) and (d) we present the corresponding phonon 
spectral functions, as obtained from ED. Again we find pro- 
nounced differences between the adiabatic [Fig.|3b)] and the 
anti-adiabatic [Fig.|2jd)] cases, revealing the distinct nature of 
the Peierls instability — driven by A or g — for small and large 

LO /t. 

For ujo/t = 0.4 [Fig.|2jc)], we observe a peak at lu = 
originating from the homogeneous (Q = 0) shift of the elec- 
tronic level for A > 0, as well as a signal located at the bare 
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FIG. 2: (Color online) Left panel [(a),(c)]: Spectral density for photoemission [A~ (k,uf); red dashed lines] and inverse photoemission 
[A + (k, uj); black solid lines] from CPT in the metallic LL phase. Energies are measured relative to the Fermi level Ef. Right panel [(b), (d)]: 
Corresponding exact phonon spectral function B(Q,uj) for the allowed wave-numbers Q of iV-site clusters. Electron and phonon spectra have 
been calculated in the adiabatic [attractive LL; (a),(b)] and anti-adiabatic [repulsive LL; (c),(d)] weak-to-intermediate EP coupling regimes. 
Dashed vertical lines in (b), (d) indicate the bare phonon frequency. 



phonon energy u>q. More importantly, the phonon excitations 
gradually soften near the zone boundary (Q — ir) already in 
the LL phase, a behavior characteristic of the Peierls transition 
traditionally regarded as a displacive phase transition. 

For large c^o A [Fig-EJd)], we observe two main absorption 
features. The first excitation near uiq becomes even harder 
with increasing wave number, in contrast to the phonon soft- 
ening in the adiabatic regime. The second branch emerging 
from u> — (Q — 0) resembles the two-spinon continuum 11 
of the XXZ model (onto which the spinless Holstein model 
may be mapped in the anti-adiabatic strong-coupling limit 2 ), 
i.e., it can be traced back to phonon-signatures of the corre- 
sponding electronic excitations. 

We now turn to the insulating state of the HMSF (Fig.|3jl. 
Due to the Peierls distortion, and the formation of long-range 
CDW order for N — * oo, a gap opens at Ep. In the adiabatic 



regime, the transition is from a LL to a traditional Peierls band 
insulator, whereas in the anti-adiabatic regime, the charge car- 
riers first undergo a cross-over to small-polarons (polaronic 
metal), which then (upon increasing the coupling further) or- 
der to form a polaronic superlattice. 

From Fig. [3ja) we see that the gap in the single-particle 
spectrum is still quite small. This is partly due to the fact that 
CPT in the form used here does not fully take into account 
long-range order. At fcp, the strongest signatures lie just above 
and below Ep, with the incoherent part again being very small. 
Away from kp, the incoherent band is significantly broadened 
(again almost oc g 2 at k — ir), and reflects the Poisson distri- 
bution of the phonons in the ground state. Nevertheless, we 
can still detect the dispersion of the split electronic band. 

In contrast to the adiabatic case of Fig.[3ja), Fig. 0c) shows 
a clear gap for all k. The polaronic charge carriers — in addi- 



4 



3 



3 

+ < 



(a) co /t = 0.4, g =5 [N =6] 



>JvV\,\ A K • 

.'.'V.A/iAii'; ■ 

-"I >'l .'J, < 





(CB-Ep)/t 



3 

< h 



3 



(c) 0) n /t = 4, g=5 [N =4] 



-30 



-20 -10 10 
(co-EJ/t 



3 






1 

(b) co Q / 1 = 0.4, 

/ 


2 ' 

g =5 [N = 6] 




















,1 


i 



0.2 



0.4 
CO/t 



0.6 



3 4 
(0/ t 



0.8 



/ 


' 1 
(d) (0 /t = 4, 

u 


2 

g =5 [N = 4] 






I 




1 

/I 
y i 




, 1,1,1,1,1,1,1,1, 



FIG. 3: (Color online) As in Fig. [2] but for parameters in the insulating CDW phase, and for the case of a band insulator [(a),(b)] and a 
polaronic superlattice [(c), (d)]. 



tion to forming a superlattice — are quasi-localized, leading to 
dispersionless excitations near energies Ep ± nuo. 

The phonon spectra for the insulating phase are shown in 
Figs. |3b),(d). In the adiabatic case of Fig. |5Jb), the zone 
boundary phonon has almost completely softened. In the limit 
TV — > oo, we expect a perfect degeneracy of Q — 0, ir for 
g > g c , and a macroscopic population of the soft phonon 
mode. For ojo/t — 4 [Fig.|3jd)], there exist two completely 
fiat signatures at uj = and ui = ujq [cf. Fig.01d)], the former 
having very large spectral weight due to the large "phonon 
content" of the polaron band. 

Similar to the results for A(k,u>), it is desirable to obtain 
results for the phonon spectrum of the infinite system. How- 
ever, in deriving a CPT-like approach to the phonon Green 
function D(q, u>), it turns out that the first-order correction in 
t vanishes identically, since the electron number per cluster 
is a conserved quantity. Therefore, CPT results in a simple 
Fourier transformation of the cluster Green function (again 



obtained by the KPM), 

1 n c 

D(q,u) = — J £Dlf(u)e- i ^-rj). (8) 
c «.i=i 

Notice that this result still becomes exact for the free Green 
function Do{q,uS), and in the limit of strong coupling and 
large phonon frequency, when t is a small parameter. As the 
calculation of the second-order contribution is extremely de- 
manding we have restricted ourselves to this expansion to il- 
lustrate the ^-dependence of B(q, uo) in Fig.|4] 

The density plots in Fig. [4] clearly summarize the differ- 
ences between the adiabatic and anti-adiabatic regimes, and 
between the LL and CDW phases of the HMSF In Fig. Ufa) 
we see the renormalized phonon dispersion uj(q), which soft- 
ens with increasing EP coupling, leading to a degeneracy of 
excitations at Q = 0, ir at g c [see panel (b)]. The strong zero- 
energy absorption feature at 7r/2 is an artifact of the small 
cluster size and the open boundary conditions used in the CPT 
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FIG. 4: (Color online) Density plots of the (CPT) phonon spectral 
function B(q,ui), where panels (a), (b), (c) and (d) correspond to 
the attractive LL (wo/i = 0.4, g 2 = 2.5), band insulator (cjo/t = 
0.4, g 2 = 5), repulsive LL (u>o/t = 4, g 2 — 0.5), and polaron 
superlattice regimes (luo /t — 4, g 2 = 5), respectively. 



scheme [cf. Fig. 0b)]. Above the Peierls transition we find — 
in agreement with recent Monte Carlo simulations 5 — that the 
soft Q — tt phonon mode splits into two branches with the 
upper one hardening as the EP coupling increases further. 
Quite differently, in the anti-adiabatic case, we observe 



two phonon signatures for all g > 0. In the LL phase 
the bare phonon mode hardens, whereas a second mode be- 
comes strongly over-damped near Q = tt [Fig. EJc)]. Finally, 
Fig- Eld) reveals a dispersionless signal at uj = u>q, as well as 
the flat polaron band at u> « for the polaronic CDW state. 
Thus, with increasing phonon frequency, we find a cross-over 
from a soft-mode (displacive) to a central-peak-like (order- 
disorder-type) phase transition, similar to the analysis of the 
spin-Peierls transition motivated by CuGe0 3 £ 



V. SUMMARY 

By presenting highly-reliable numerical results for the 
electron- and phonon spectral function, we have identified 
four physically distinct regions in the phase diagram of the 
one-dimensional spinless Holstein model, characterized by 
important generic features such as attractive or repulsive Lut- 
tinger liquid behavior and phonon softening, which are highly 
relevant for low-dimensional condensed matter systems. 
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